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Abstract 

The aim of this paper is to find out how would possible space 
non-commutativity (NC) alter the QM solution of the Coulomb prob- 
lem. The NC parameter A is to be regarded as a measure of the 
non-commutativity - setting A = means a return to the standard 
quantum mechanics. As the very first step a rotationaly invariant 
NC space R|, an analog of the Coulomb problem configuration space 
Rq = R 3 \ {0}, is introduced. R 3 is generated by NC coordinates 
realized as operators acting in an auxiliary (Fock) space T . The prop- 
erly weighted Hilbert-Schmidt operators in T form %\, an NC analog 
of the Hilbert space of the wave functions. We will refer to them as 
"wave functions" also in the NC case. The definition of an NC ana- 
log of the hamiltonian as a hermitian operator in H.\ is one of the 
key parts of this paper. The resulting problem is exactly solvable. 
The full solution is provided, including formulas for the bound states 
for E < and low-energy scattering for E > (both containig NC 
corrections analytic in A) and also formulas for high-energy scattering 
and unexpected bound states at ultra-high energy (both containing 
NC corrections singular in A). All the NC contributions to the known 
QM solutions either vanish or disappear in the limit A — > 0. 
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1 Introduction 



Basic ideas of non-commutative geometry have been developed in [T] and, in a 
form of matrix geometry, in [2] . The main applications have been considered 

• in the area of quantum quantum field theory in order to understand, or 
even to remove, UV singularities, and eventually, 

• to formulate a proper base for the quantum gravity. 

The analysis performed in [3] led to the conclusion that quantum vac- 
uum fluctuations and Einstein gravity could create (micro) black holes which 
prevent localization of space-time points. Mathematically this requires non- 
commutative (NC) coordinates x^ in space-time satisfying specific uncer- 
tainty relations. The simplest set of NC coordinates x M should satisfy Heisen- 
berg-Moyal commutation relations in an auxiliary Hilbert space 

[a^ar"] = iO^, fi, v = 0, 1, 2, 3 , (1) 

where 9^ u are given numerical constants that specify the non-commutativity 
of the space-time in question. 

Later in jl] it was shown that field theories in NC spaces with (JTJ can 
emerge as effective low energy limits of string theories. These results sup- 
ported a vivid development of non- commutative QFT. However, such models 
contain various unpleasant and unwanted features. The divergences are not 
removed, on the contrary, UV-IR mixing appears, [5 J . The Lorentz invariance 
is broken down to 5*0(2) x SO(l, 1), but even this is sufficient to prove the 
classical CPT and spin-statisics theorems, j6]. This was not accidental and 
led to the twisted Poincare reinterpretation of NC space-time symmetries, 
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However, it could be interesting to reverse the approach. Not to use the 
NC geometry to improve the foundations of QFT, what is a very compli- 
cated task, but to test the effect of non-commutativity of the space on the 
deformation of the well-defined quantum mechanics (QM): 

• Various QM systems have been investigated in 3D space with Heisen- 
berg-Moyal commutation relations [a^a^] = i9 1 ^ , i,j = 1,2,3, e.g. harmonic 
oscillator, Aharonov-Bohm effect, Coulomb problem, see [8], [9]. However, in 
such 3D NC space the rotational symmetry is violated and there are systems, 
such as iJ-atom, that are tightly related to the rotational symmetry. 

• The rotational symmetry survives in 2D Heisenberg-Moyal space with 
NC coordinates £1,22 satisfying the commutation relations [a; 1,0:2] = iO in 
an auxiliary Hilbert space. In [10] a planar spherical well was described in 
detail: 

(i) First, the Hilbert space H. of operator wave functions ip = tp(xi,X2) 
was defined; 

(ii) Further, the Hamiltonian was defined as an operator acting in 7-L. It 
was nice to see how the persisted rotational symmetry helps to solve exactly 
the problem in question. 

The presented list of references is incomplete and we apologize for that. 
We restricted ourselves to those which initiated progress or are close to our 
approach. 

Our aim is to extend this scheme to the QM problems with rotationally 
symmetric potentials V(r) in the configuration space i?g = R 3 \ {0}. We 
restrict ourselves to the Coulomb potential which, in the usual (commutative) 

3 



setting, is a solution of the Laplace equation vanishing at infinity: 

AV(r) = V(r) = ~\- ( 2 ) 

For Coulomb problem , in a Gaussian system of units, q is a square of electric 
charge q = ±e 2 : q > or q < corresponding to the Coulomb attraction 
or repulsion respectively. We are dealing with Schrodinger equation 

h 2 a 

- — A^(x) - -V>(x) = Eip(x), r=|x|>0 (3) 
2m r 

in the Hilbert space T-L specified by the norm 

Ml = J rf 3 x|^(x)| 2 . (4) 

Expressing the wave function as 

^(x) = Rj(r)H jm (x), HjmW ~ T*Y jm (ti,ip), (5) 

and putting a = mq/h 2 and k 2 = 2mE/h 2 , we obtain the radial Schrodinger 
equation: 

r B!-(r) + (2j + 2)i2j(r) + 2ai? i (r) = - fc 2 r ^(r) . (6) 

The parameter a is simply related to the H-atom Bohr radius a® = h 2 / (me 2 ) = 
\at\~ 1 . The solutions of (|6]) are given in terms of solutions of the confluent 
hypergeometric equation (see, e.g. [TT]). 

(i) For negative energies one obtains bound states with discrete energy 
eigenvalues: 
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Tfl (■ 

En = -^ 2 , n=j + l, 3 + 2,.... (7) 
4 



(ii) For positive E 
partial wave S-matrix 



— k 2 

2m 



> one obtains scattering states with j-ih 



Eg + 1 ~ I ) 

r (j + i + f ) ' 



(8) 



This paper is organized as follows. We define the NC QM Coulomb problem 
in Section 2 : (i) We define the rotationally invariant NC configuration space 
R| - the NC analog of the Coulomb problem configuration space Rq and 
the Hilbert space %\ of wave functions in R|, then, (ii) we introduce an 
important hermitian operator - the NC Coulomb problem Hamiltonian. In 
Section 3 we derive the NC analog of the radial Schrodinger equation, and 
we solve it for all energies. For energy E < we recover NC bound states 
regular in our non-commutativity parameter A, while for E > one obtains 
two sectors - low energy scattering regular in A, and for ultrahigh energies 
there are solutions singular in A which disappear in the commutative limit 
(this extends our partial results for E < in p2]). The last Section 4 contains 
discussion and conclusions. 



The name "non-commutative quantum mechanics" may seem to consist of 
more words than neccesary, since one cannot think of a quantum theory lack- 
ing certain non-vanishing commutators - the underlying uncertainty cannot 
be detached from the theory. The way from QM to NCQM can be roughly 
described as an analogy of the transition from the classical theory to QM. 
In the latter one the Heisenberg principle results into fuzziness of the phase 



2 The non- commutative space R x 
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space. In NCQM even the notion of a single point in the configuration space 
loses relevance. This fact is reflected in the non-vanishing commutator of the 
coordinates involved. Since we are about to deal with the Coulomb problem, 
the commutation relations have to preserve the rotational symmetry. 

In this section we define the non-commutative (NC) space R|, possessing 
full rotational invariance, as a sequence of fuzzy spheres introduced, in various 
contexts, in [13]. Different fuzzy spheres are related in such a way that at 
large distances we recover space Rq with the usual flat geometry. A similar 
construction of a 3D NC space, as a sequence of fuzzy spheres, was proposed 
in [13] . However, various fuzzy spheres are related to each other in a different 
manner (not leading to the flat space Rq at large distances). 

The non- commutative configuration space 

We realize the NC coordinates in R| in terms of 2 pairs of boson annihilation 
and creation operators a a , aj,, a = 1, 2, satisfying the following commutation 
relations, see [T5] : 

[a a , ap] = <W> [ a a, a p\ = l a L a l] =0. (9) 

They act in an auxiliary Fock space J 7 spanned by normalized vectors 

(a\) ni (4) n2 . . 
|m,n 2 ) = 1 1, r l2L . (10) 

Here 1 0) = |0, 0) denotes the normalized vacuum state: a\ |0) = 02)0) = 0. 

We shall use the notation T n = {\n\, n 2 ) | n\ + n 2 = n}. 

The noncommutative coordinates Xj, j = 1,2,3, in the space R^ are given 

as 

Xj = Xa+cTja = \cr J a[S alap, j = 1,2,3, (11) 
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where A is a universal length parameter and aj are Pauli matrices. The 
operator that approximates the NC analog of the Euclidean distance from 
the origin is r = A (N + 1), N — a^ a a a . The coordinates Xj and r satisfy 
rotationally invariant relations: 

[xi,Xj] = 2i\e ijk x k , [x u r] = 0, r 2 - x] = A 2 . (12) 

We will provide a strong argument supporting the exceptional role of r later. 

Hilbert space Hx of NC wave functions 

Let us consider the linear space of normal ordered analytic functions con- 
taining the same number of creation and annihilation operators: 

* = Yl C m m n 2 {a{) mi (4) m2 («l) ni , (13) 

where the summation is finite over nonnegative integers satisfying m 1 + m 2 = 
m + n%. Here %\ is our denotation of the Hilbert space of functions ( TT3]) 
which possess finite weighted Hilbert-Schmidt norm 

||*|| 2 = 4vrA 3 Tr[(iV + l)^ t *] = 4vr A 2 Tr[r ^ . (14) 

The rotationally invariant weight w(r) = 4tt A 2 r is determined by the re- 
quirement that a ball in R| with radius r should possess a standard volume 
in the limit r — > oo. The projector P n on the subspace Fq@ . . . (BJ- n , corre- 
sponds to the characteristic functions of a ball with the radius r = X(N + 1). 
Therefore, the volume of the ball in question is 

n . , 

V r = 4vrA 3 Tr[(iV+l)PJ = 4vr A 3 V(Jfe + l) 2 = — r 3 + o(-) . (15) 

Thus, the chosen weight w(r) = 4tt X 2 r has the desired property. 
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Orbital momentum in H\ 

In 7-L\ we define orbital momentum operators, the generators of rotations Lj, 
j = 1, 2, 3, as follows 

LJ* = ^[a + a ja ,% j = 1,2,3. (16) 

They are hermitian (self- adjoint) operators in 1-L\ and obey the standard 
commutation relations 

[L^LjYSl = (LiLj - LjL^ = ie ijk L k ty . (17) 

The standard eigenf unctions ^j m , j = 0, 1, 2, . . . , , m = —j, . . . , +j, sa- 
tisfying 

Ll^ jm = j(j + l)V jm , L 3 ^ jm = m^ jm , (18) 



are given by the formula 

mi! m 2 ! 3KtiJ nil n 2 ! 



v = E (**„, b) <W!, (19) 

J 111 .1 lll,,t II ,' II.J 

(jm) 

where g = Xa* a a a = XN. The summation goes over all nonnegative integers 
satisfying mi + m 2 = ni + n 2 = j, m\ — m 2 — n\ + n 2 = 2m. Thus \^- m = 
when restricted to the subspaces T n with n < j (what corresponds to the 
fact that in the standard QM the first j — 1 derivatives of ^ j m vanish at the 
origin). For any fixed Rj(g) equation f|T9|) defines a representation space for 
a unitary irreducible representation with spin j. 

The radial part and normalization in 7-L\ 

The two wave functions ^ j m and ^j' m ', with (j,m) ^ (j',m') and arbitrary 
factors Rj(g) and Rj'(g), are orthogonal in 7-L\. Thus, when evaluating the 
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norm of ^j m , it is sufficient to calculate j | "^1/^^ 1 1 2 = H^jH 2 (this equality 
follows from the rotational invariance of the norm in question): 



oo n 



\\V jm \\ 2 = 4ttA 3 + {k,n-k\{n + l)^]^ j:j \k,n-k) . (20) 

n=j k=0 

We benefit from the fact that tyjj has a simple form 

*ii = (a\y R^e) (-a 2 y . (21) 

The matrix element we need to calculate is 

(k, n — k\ (a\y Rj(g) a{ (a\y Rj{g) a\ \k,n — k) 

(k + j)l(n-k)\ 
k\ (n — j — k)\ 

where 



|72i(n-j)r, (22) 



R.{n) = (k,n- k\Rj(Q)\k,n- k) (23) 

(the expression on the r.h.s. is k - independent). Inserting (l2"Tj) . (I22p into 
(|2"U|) and using the identity (see [T7]) 



£ 

fc=0 

we obtain 

oo 



k + j\ fn — k\ fn + j + 1 



3 J \ 3 J V 2j + l 



ll9 4ttA 3+2j ' , /n + 2?' + l\ . . 



(jiy ^ v 2i+i 

This expression represents, up to an eventual normalization, the square of a 
norm of the radial part of the wave function. 

The NC analog of Laplace operator in 7-L\ 
We postulate the NC analog of the usual Laplace operator in the form: 

Aa^ = - y r [al [a a , n = - * + [al [a a , *]] . (25) 



This choice is motivated by the following facts: (i) A double commutator is 
an analog of a second order differential operator, (ii) the factor r _1 guarantees 
that the operator Aa is hermitian (self-adjoint) in TL\, and finally, (iii) the 
factors A" 1 , or A~ 2 respectively, guarantee the correct physical dimension of 
A x and its non-trivial commutative limit. 

Calculating the action of (j2"5]l on ^ j m given in ( TT91 we can check whether 
the postulate (125]) is a reasonable choice. 

The operator Rj{g) in (1191) can be represented as a normal ordered expansion 
of an analytic function TZj(g) : 

We) =:K*G?):=£c> :g k := £4^—^— (26) 

k k 

The last equality follows from the equation 

n\ 

■ N : \ni,n 2 ) = t -ttt \ni, n^), n = n 1 + n 2 (27) 

(n — ky. 

(which can be proved by induction in k). Since : N k : \nx,n2) = for 
k > ni + ri2, the summation in fl26|) is effectively restricted to k < n on any 
subspace T n . 

The following formula follows from commutation relations (186]) ; the proof 
is given in Appendix A: 

[4, [a a , %\ = X> ^ {alr \ {al r X 

(jm) 

n ni ( —n ) U2 

x : \-qW\q) - 2{j + l)K>{e)\ ■ 1 \ 2 ■ (28) 

nil n 2 \ 

Here TZ'(g) denotes the usual derivative defined as: 

oo oo 
fc=0 fc=l 
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and 1Z"(g) is denned as the derivative of TZ'(g). Thus, the prime corresponds 
exactly to the usual derivative d e . In the commutative limit A — > operator 
g formally reduces to the usual radial r variable in R 3 , and we see that A a 
just reduces to the standard Laplace operator in R 3 . 

The potential term in 7-L\ 
The operator V corresponding to a central potential in QM is defined simply 
as the multiplication of the NC wave function by V(r): 

(VV)(r) = V(r)V = tfV(r). (30) 

Since any term of \& G contains the same number of creation and anni- 
hilation operators (any commutator of such a term with r is zero), the left 
and right multiplications by V(r) are equal. 

In the commutative case the Coulomb potential is a radial solution of 
the equation (j2J) vanishing at infinity. Due to our choice of the NC Laplace 
operator A\ the NC analog of this equation is 

A x V(r) = & [at, [a a ,V(N)}} = 0. 

Last equation can be rewritten as a simple recurrent relation 

(N + 2)V(N + 1) - (JV+1) V(N) = (N + 1)V(N) - NV(N-l) . (31) 

Putting (M + l)V(M) - MV(M - 1) = g and V(0) = g - f, and 
summing up the first equation over M = 1, ... N, we obtain the general 
solution: 

V W = -xWTT) + » = -* + *• (32) 
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where q and go are arbitrary constants (A is introduced for the future conve- 
nience). Thus the NC analog of the Coulomb potential vanishing at infinity 
is given by 

(Vtf)(r) = (33) 
r 

We see that the ~ = X {n+i) dependence of the NC Coulomb potential is 
inevitable. 



3 The Coulomb problem in NC QM 

Based on f )25|) and fl33|) we postulate the NC analog of the Schrodinger equa- 
tion with the Coulomb potential in R| as 

[al [a a , n--^ = E^ & \ [al [a a , *]] - 2a $ = k 2 , (34) 



In Appendix A these two equations are proved: 

[at, [a a , * jm }} = J2 ■■■ ■ i-e Xn " ~ 2(j + l)A^] : ... , 

r * im = ■ ■ ■ : KQ + X J + A) TZj + XgK'j] : ... , (35) 

U m ) 

where TZj = TZj(g) and similarly for derivatives. The dots on the left and 
right in ( 135]) denote the products in \1/ containing respectively creation and 
annihilation operators together with the factor A- 7 , that represent the angular 
dependence of \l/ and remain untouched as the operators in question are 
rotation invariant. Inserting (135]) into ( I34|) we obtain the NC analog of radial 
Schrodinger equation: 

: qTZ] + [k 2 \g + 2j + 2] K'j + [k 2 Q + k 2 X(j + 1) + 2a] TZj : = . (36) 
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We claim ( 13 6 p to be an NCQM analog of the usual radial Schrodinger 
equation ([6]) known from QM. There definitely is a resemblance; in the limit 
A —7- the terms in fl36|) proportional to A representing the NC corrections 
disappear. Considering the same limit we see that the presence of the colon 
marks denoting the normal ordering should not worry us either; recall that 
for zero A it makes no difference whatsoever whether we care for the ordering 
or not. (Normal and usual powers coincide for zero A.) This is a good news 
to start with, leaving us, however, with the task to solve ( 13 6 p for nonzero 
A, which means that both the extra terms proportional to A and the normal 
ordering are to be taken at a face value. If it was not for the ordering issues, 
the solution would be quite straightforward - the extra terms would mean 
just adding some more work needed to complete the calculation, but it is 
known how to solve the problems of this kind. In fact this is precisely what 
we are going to do: We associate the following ordinary differential equation 
to the mentioned operator radial Schrodinger equation 

gK'j + [k 2 Xg + 2j + 2] IZ'j + [k 2 g + k 2 X(j + 1) + 2a] Kj = 0, (37) 

with q being real variable, and we will solve this one. But how come we 
expect this step to be of any use to us, when we actually do have to care 
about the ordering ? One should notice the following: whatever appears in 
the equations, it can be expressed in terms of powers in g ; it is just that 
they are normal powers in one case and the usual powers in the other. Next, 
we have some operators in both equations; normal derivatives in one case 
and the usual ones in the other. The key information is, that the derivative 
defined in ( 129]) acts on the normal powers just like a carbon copy of usual 
derivative, see ( 1S5|) . 
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Now bearing this in mind, we expect R = : 1Z : , the solution of (135)) . to be of 
the same form as 71, the solution of (137)) . except for the nature of the powers 
involved. So a brief summary goes like this: The solution of ( 137)) with all the 
usual powers replaced by the normal ones is the solution of ( 1361) . However, 
the form of the solution is not the best one yet. We have already mentioned 
the relation between the equation given by QM (j5J) and ( 1361) . the one supplied 
by NCQM. Of course we would like to compare the corresponding solution 
as well, but this is rather a difficult task as long as we have normal powers in 
the first and the usual ones in the latter one. Fortunately we have a formula 
relating : g n : and g n , namely 

Nl iV' 

n. = X n. N n ._ . -n . = A ~n . jyn . = X ~n /g g x 

(N - n)\ (N - n)\ 

All we need is to rewrite : 7Z : using those relations. Then the above men- 
tioned comparison of QM and NCQM will be obtained. 

Perhaps is has been made clear enough what is to be done, so let's get started 
with the solution of equations ( 137)) and (1361) . Some mathematic theory is to 
be studied here. We have got a second order differential equation of the form 

(a x + b ) y"(x) + (a x x + b x ) y'(x) + (a 2 x + b 2 ) y(x) = (39) 

Depending on whether the quantity D 2 = a\ — Aaoa 2 is zero or not, the 
solutions of ( 139)) are given in terms of Bessel or confluent hypergeometric 
functions respectively. We will restrict ourselves to the case ao = 1, b = 
at the price of some generality loss -but generality is not what we are after 
in the first place. 
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(i) a — 1 , b = 0, D 2 = a\ - 4a a 2 7^ 

In this case, the solution of (13T?j) is of the form 

= e £ ^ X Z{a,c,-Dx), (40) 

where a = ( D 2 ai ^i + ^2) , c = 61, and Z(a,c,x) is any solution of the 
confluent hypergeometric equation 

xZ' + (c-x)0 - aZ = 0. (41) 

The solution of (|4"T]) regular at the origin will be the most important to us. 
It is often referred to as the confluent hypergeometric function: 

4,(a;c;S) = £) (42) 

m=0 ^ m m ' 

where (a) m denotes the so-called Pochhammer symbol: (a)o = 1 and 

r(a + m) 

(a) m = a(a + 1) ... (a + m-1) = — =^-r — , m = 0,1,2,... . 

T(a) 

The fundamental system of solutions of (|41[) consists of the following func- 
tions: 

if) (a, c; x) , and e x if)(c — a, c; — x) , (43) 
with if>(a,c;x) possessing this asymptotic expansion for x — > 00 : 

if)(a,c;x) = V(-ir (a)m(a ~, C+1)m ^ m . (44) 

Every solution (e.g. also (|42p . the one regular at the origin) can be expressed 
as a suitable linear combination of (j4"3"|) . This possibility is useful when treat- 
ing the scattering processes . 
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Note that D is fixed up to the sign, since the coefficients in the equation 
determine the value of D 2 only. In fact it does not matter which possibility 
is preferred in pUj) . Replacing D with — D makes no difference because of 
the Kummer identity which reads 

e~ x/2 4>(a, c; x) = e x/2 0(a, c; -x) . (45) 



(ii) a = 1 , b = 0, D 2 = of — 4a a 2 = 

This time the solution of fl39|) has the following form: 

y (x) = e-^'x^C^ (y{-2a l b l +Ab 2 )x) (46) 

C v (x) is any solution of the Bessel equation 

x 2 C'l + xC' u + (x 2 - v 2 )C v = (47) 

We require regularity in the origin, what leads us to the Bessel function: 

£W*) = J.(*) = E(-i)-(|) m ,r (m + , + i) < 48 ' 



m=0 



After this general mathematic interlude we return to our equation (|37|) . For 
the sake of brevity it is suitable to introduce a new parameter r\ defined as 



kX V2mEX /A , 

. = T = ^— (49) 



According to (1371) the following holds: 

J 4i? 2 

(50) 



a = 1 ai = \k 2 = ^ a 2 = A; 2 = ^ 



6 = &! = 2(j + 1) 6 2 = A/e 2 (j + 1) + a = |(j + 1) + 2a 
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D = ±(\ 2 k 4 - Ak 2 ) 1 ' 2 = ± T vV (V 2 ~ 1) 

A 

For D ^ 0, or equivalently r) ^ 0, 77 7^ 1 the solution of (l3T)j) is 

Rj± = : 72.j ± : 
= : exp 



_ | _ 2??^/r ? 2~l _ 2r^_ 



X 



X H J + 1± ^7?=T' 2J + 2 ' T ^) : 



2^2-1 

...see Appendix... 
1 ± 2r]^r] 2 -1-2t] 2 



(51) 



1 W 



xF [j + 1 ± 



Q:A 



-iV, 2j + 2, 



4r?Vr? 2 -l 



'l±2r? A /r? 2 -l-2r; s 



F(a, 6; c; z) is the usual hypergeometric function: 

{a) m {b) m z r ' 



F(a,b; c; z) = 



m=0 



(c) m m! ' 



It is one of the solutions of the hypergeometric equation 

z{\ - z)u"{z) + [c - (a + b + l)z]u\z) - abu(z) = 0. 



(52) 



(53) 



The calculations needed to get rid of the normal ordering in the above equa- 
tion are briefly sketched in Appendix. The ± signs that emerged as a lower 
index in R± spring from the two possible choices of the sign of D. We have 
mentioned that the choice of sign is completely arbitrary due to the Kummer 
identity which holds for the confluent hypergeomatric function. This fact 
survives the process of rewriting the normal powers in terms of the usual 
ones and is reflected in the so-called Euler identity for the hypergeometric 
functions: 



F(a, b, c; x) = (1 — x) b F(c — a, b, c; x(x — 1) x ) . 



(54) 
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If rj = 0, the solution of fl36|) is 



= : Q-i-W J_ 2i _ x (V8^) 
= ...see Appendix... 



(55) 



■ Q^^ft-N, 2j + 2, 2aX) 



And finally for 77 = 1 



i?^ - : TZj : 

= : e WA^-i/2j_ 2 ._ l(v ^). 

= ...see Appendix... 



(56) 



To sum up, the solution of (1361) typically consists of an exponential factor 
(or its NCQM analog, see Appendix for further explanation) multiplied by 
certain power series. Depending on the energy, the exponent may be real or 
imaginary. We are supposed to investigate the following cases: 

(i) 77 G iR or E < 0; 

(ii) V G (1, +00) or £>S; 
(hi) 77 G (0,1) O r0<£<f|; 
(iv) 77 = or E = 0; 

In (i) and (ii) the arguments in the exponentials are real - bound states 
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may occur for certain energy values. The (iii) leads to the scattering, the 
exponential factor involves an imaginary part. The last two, (iv) and (v), 
are the border cases separating the above mentioned intervals. 



3.1 Bound states 

At the very beginning it is suitable to remind us of the QM version, which 
predicts bound states to occur under the condition that the potential is 
attractive (a > 0) and the energy has some specific negative values mentioned 
in (ED, i.e. 

4 

me 

E n = 77^2 2 ' n = 3 + !> 3 + 2 ' ■ ■ ■ ■ 

Now to NCQM. In this section, the R + form of the solution of (155]) will be 
suitable for our purposes. The reason is, that in R + the absolute value of 
the factor multiplying the hypergeometric function is less than 1 for every 
N. Consequently, when looking for bound states, is is sufficient to check 
whether the power series terminates. In case of hypergeometric function this 
happens if the first argument is a negative integer. This leads to discrete 
energy values. 

3.1.1 Bound states for E < , rj — i \r}\ 
The equation floT]) in this case reads 



Rj 



l-2Mv/M 2 + l + 2| 



A' 



x 



xF \j + 1 . -N; 2j + 2; J j== 

(57) 
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Square integrable bound state solutions can be easily seen, since the hyperge- 
ometric function has to reduce to a polynomial. This happens only if a > 0, 
i.e. in the Coulomb attractive case, provided that 

aX , , 

a > and 7 + 1 . = —n . (58) 

2|77|^tF + T 

This gives the bound state energies (with Planck constant H and mass m 
explicitly introduced): 

E[ n = Jgl 2 = f 1 - + («A/n) 2 ) . (59) 

Taking the limit A — > we recover the QM result. We see that the NC 
corrections are governed by (Xa) 2 = (X/a ) 2 - the square of the ratio of two 
dimension-full parameters in the model - the noncommutativity parameter A 
and the Bohr radius a Q . 

3.1.2 Bound states for E > 2/X 2 , rj = \rj\ > 1 
The equation (J5"T]) in this case is 



Rj 



l + 2|77|v/|77| 2 -l-2| 



JV 



xF| J + l+ , , a " .,-iV;2 J + 2; 4| ^V^ 



2\ V \^\ V \ 2 - l' l + 2\r]\^\T]\ 2 -l-2\T]\ 



2 



(60) 

Since the absolute value of the prefactor preceding the hypergeometric 
function is less than 1, the whole solution can have a finite norm. There 
is a possibility for the hypergeometric function to terminate, since the first 
argument becomes a negative integer for certain energy values, under the 
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condition that a < (so the potential has to be repulsive this time). 

ct 

a < and j + 1 H . = — n . (61) 

In this case the bound state energies read: 

, I _ _ _ x 2h 2 



E 



n 



(l + y/l + (aX/n) 2 ) = ^-E{ n . (62) 



An " mX 2 

These are very unexpected solutions, E^ n being a mirror of E{ n with respect 
to the "critical energy" E crit = 2h 2 /(mX 2 ). However, they disappear from 
the Hilbert space 7i\ in the commutative limit A — > 0. 

For energies E = E\^ n the solution (1601) has the same finite norm as (1571) for 
E = E{ n . 

3.2 Scattering E G (0, 2/A 2 ), rj G (0, 1) 

In this section we will deal with the NCQM version of Coulomb scattering. 
We will briefly sum up the QM results before handling our NCQM case. 
The solution of the radial Schrodinger equation for a particle in a potential 
V(r) = —a/r with the angular momentum j and energy E > 0, regular in 
r — > is given in terms of the confluent hypergeometric function (see [TT]): 



R 



QM 



e ikr + ij, 2j + 2, -2iAx) , k = V2E > . (63) 



We have refrained from writing down m/h 2 explicitly here. This will simplify 
the formulas and will not do any harm, since the full form can be restored 
anytime. 

The solution is real and for r — > oo it can be written as a sum of two complex 
conjugated parts corresponding to an in- and out- going spherical wave. In 

21 



the following formula a real term common for both parts is left out, having 
no influence on the S-matrix. 



J)QM . ifcr+if ln(2fcr) , 6 _~ifcr-i§ ln(2fcr) / fi4 \ 

r(j + i + if) e + ru + i-q) • lb4j 

The S"-matrix for the j-th partial wave is defined as the ratio of the r- 
independent factors multiplying the exponentials with the kinematical factor 
left out. 

Now let us have a look on the Coulomb scattering in NCQM, consid- 
ering the j-th partial wave and the energy E G (0, 2/A 2 ). For the future 
convenience it is suitable to introduce 



p = ^2E{l- l -\*E) . (66) 

The solution regular in the origin is given in terms of the hypergeometric 
function (152]) : 

i?£y = - — ^ F j + l-i- -iV ) 2j + 2;2iAp^— — 67 
\p — i\E) \ p p + iXEJ 

We choose the positive square root (jB"6"|) for £ G (0, 2/A 2 ). The formula 
(l66|l represents a conformal map from the upper complex E'-plane on a right 
complex p-plane with a branch cut p G (0, 1/A). The radial dependence 
of Rj is present in the hermitian operator iV: r = g + A, g = XN. In 
analogy with (J2H ) we will rewrite also the NC solution as a sum of two 
terms corresponding to the in- and out- going spherical wave. Again, leaving 
out the common hermitian factor which is irrelevant regarding the S'-matrix, 
we can write (see Appendix): 
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j+1 oar/p 



Rej ~ r(; + i + <f) Vp-*A£ 



p + iXE 



N+j+l—i 



* T(2j + 2)T(N+1) 
T(N + 2+j - ia/p) 



(2Ap) 



. .a „ T „ . .a — z fp + iXE 
x F j + l-i- -j-i- iV + 2 + j-i-; — - — — 
p p p IXp \p — iXE 



+ 



r(j + i-if) Vp + ^ae 



P-^r +j+1+i " r(2 J + 2)r(iv + i) 

J r(iv + 2 + j+Wp) l2Apj 

p — iA-E 1 



x F(j + l + i-,-j + i-,N + 2 + j + i-,— ( F )). 6^ 

\ p p p 2Ap \p + zA-fc J J 

To enable better comparision with f lM|) let us rewrite also ([55]) as a sum of 
two complex conjugated parts. Some sort of sketch of the calculation leading 
to it is to be found in the Appendix C: 

r(2j + 2) e -^ ln(2pr) 



R 



Ej 



j+l — j-l -aiv/2p 



x exp 
x exp 



— (r/A + j + ia/p) In 



r(j + 1 - ia/p) (2pr)J +1 
p + iXE 



p — iXE 



n=l 



B n+1 (j + l + ia/p) - B n+1 (0) 
n(n + 1) 



a .a r . „ .a 1 p — iXE 

x F[j + l + i-, -j+l-, + l + -— - F 

p p X p 2%Xp p + iXE 



+ 



x exp — (r/A + j — ia/p) In 

oo 

x oxi > -$^(V 



r(2j + 2) 



if ln(2pr) 



T(j + 1 + ia/p) (2pry +1 
p — iXE 
p + iXE _ 

n B n+l (j + l-ia/p) - B„+i(0) 



n=l 



n(n + 1) 



, a a r a 1 p + iXE 

X F [j+l-t-, -J-1-, -+J + 1-Z-; _ 

p p A p 2zAp p — zA-fe 



(69) 
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The S-matrix is the ratio of the r-independent factors : 



5? 



(E) 



r(j + 1 - »g) 
rO' + i + if) ' 



(70) 



where 



E 



1 

A 2 



1 + z v^AV - 1 



') 



(71) 



is the conformal map inverse to flBB|) . which maps the cut p right-half-plane 
into the E upper-half-plane. We take the positive square root in (I7T|) for 
p G (1/A, +oo). 

The physical- relevant values of the S'-matrix are obtained as Sj(E + ie) in 
the limit e — > + . 

The interval corresponding to the scattering E G (0,2/A 2 ) is mapped onto 
the branch cut in the p-plane as follows: 

E G (0, 1/A 2 ) i— >■ upper edge of the branch cut p G (0, 1/A) 
E G (1/A 2 , 2/A 2 ) ^ lower edge of the branch cut p G (0, 1/A) (72) 

3.3 Bound states revisited - poles of the S-matrix 

Like in the previous section, we will briefly remind the QM case: In the case 
of an attractive potential (a > 0) the S'-matrix ( )65|) has poles in the upper 
complex fc-plane for 



71 



a 
n 



n = j + 1, j + 2, 



(73) 



The wave function fl63|) is integrable for k = k n : 




(74) 



24 



It is obvious that the energy levels correspond to the poles of the S'-matrix: 

a 2 

E n = -^2, n = j + 1, j + 2, ... (75) 

In NCQM there is an analogy, the poles of the S-matrix occur in the case 
of attractive potential (a > 0) for some special values of energy below 0. 
However, poles can be found also in case of repulsive potential (a < 0) for 
particular values of energy above 2/A 2 ). 

3.3.1 Poles of the S-matrix for attractive potential 

Pn = A «>0 & E[ n = 1 (l - v/l + (Aa/n) 2 ) < 0. 

n = j + 1, j + 2, ... (76) 

In the limit A — > this coincides with the standard self-energies of the hy- 
drogen atom fl75l) . Let us denote 

n 

Then the solution (1671) is 

Rl 3 = (nifFi-n, -N, 2j + 2; -2 Kn (^)- 1 ) . (78) 

It is integrable since Q n G (0, 1) for positive k and under given conditions the 
hypergeometric function is a polynomial. The norm (j24p of R ! n j is finite and 
given in terms of a generalized hypergeometric function. We do not present 
the corresponding cumbersome formula as it is not needed for our purposes. 
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+ + Ki + 1 



K n + y/1 + K\ ~ 1 



(77) 



3.3.2 Poles of the S -matrix for repulsive potential 

(These disappear from the Hilbert space of the physical states in the limit 
A ->• ) 

P X n = ^, « < & E{ T n = ± (l + v/1 + (Xa/n)^ > 2/A 2 

n = j + 1, j + 2, . . . (79) 

Now (157)1 has the form 

Rnj = (-n I n I ) N F(-n,-N,2j + 2;2K n (Q I n I )~ 1 ) (80) 

where 

Oil /t» + y/1 + < + 1 
»n u , 2 7- 

K n - a/1 + < - 1 

The definition of K n is the same as in (1771) (note that it is negative this 
time). Since Q^ 1 = fl^ 6 (0,1) the solution ( )67j) is integrable because the 
hypergeometric function terminates like in the previous case. 



4 Comments and conclusions 

The aim of this section is to sum up what has been done, compare QM and 
NCQM results, and finally outline what could be done in the future: 
We carefully defined the NC rotationally invariant analog of the QM con- 
figuration space and the Hilbert space of operator wave functions . The 
central point of our construction was the definition of A^, the NC ana- 
log of Laplacian, supplemented by a consequent definition of the weighted 
Hilbert-Schmidt norm and a definition of the Coulomb potential satisfying 
NC Laplace equation. 
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Then we were able to introduce the NC analog of the Schrodinger equation, 
taking advantage of the spherical symmetry of the problem when separating 
the radial part. The knowledge of how are normal ordered powers of the "NC 
radial variable" correlated with the "usual" powers enabled us to solve the 
NC problem using the associated ordinary differential equation. 
Now let us see how NCQM matches the "standard" quantum mechanics. 
We are sorry to bother the reader with the QM facts which are undoubtedly 
familiar to them, but there seems to be no better way to compare the theo- 
ries. 

The quantity labeling the solutions of Schrodinger equation is energy. Some 
of the labels are excluded in the sense that they cannot be attributed to 
a physical state. In standard QM, the solutions with negative energy were 
dismissed as lacking the physical interpretation, except for those with some 
special energy values for which the wave function was normalizable. One 
would simply expect a particle with negative energy to be trapped in certain 
region, and, on the contrary, since a particle with positive energy can be 
pretty much anywhere, it is quite acceptable that Schrodinger equation does 
not provide us with normalizable solutions for a given case. 

The state of affairs seems to be a bit different in NCQM. Although the 
energy spectrum for an electron trapped in the atom is predicted in agree- 
ment with QM with small correction of order A 2 , there are two special values 
of energy, E = and E cr u = If? j (m\ 2 ) with the following feature: certain 
energy values below E = for an attractive Coulomb potential and certain 
values above E crit for repulsive potentials provide a normalizable state. 
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There is a remarkable symmetry between normalizable states corresponding 
to E ! n < and those corresponding to > E crit : 

The bound state energies are symmetric with respect to the energy \E cr a = 
h 2 /(m\ 2 ), 

EH 1 = \Ecr« (l =F y/1 + (Aa/n)=») . 

Moreover, the corresponding radial wave functions are equal up to the change 
a — > —a and the sign changing factor at each step A in the radial direction, 

= (-If 

The same symmetry can be seen for scattering states for energies: 



E 1 = ^E crit - e , p 1 = +J ^E crit - A 2 , 



E n = -E crit + £ , P H = — J -jE crit — A 2 , 



withe e (0, \E crii ). Namely, 



= (-1)" 



This relation follows directly from 



p 1 + iXE 1 _ p u + i\E u VI + A% + «Vl + A% 
pT—j^ : ~p~ir—XEii ' ' VI + A 2 £ - iVl + A 2 £ 

It would be highly desirable to see the background of those almost perfect 
reflection symmetry in the energy with respect to \E crit . 



The reader may wonder how big the NC corrections actually are. The an- 
swer is, that the parameter A is not fixed within our model. However, it 
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can be estimated by some other physical requirement. For example, one can 
postulate, as was done in early days of modern physics, that the rest energy 
mc 2 of electron is equal to the electrostatic energy of its Coulomb field. In 



this means: 



where 



8n J 



e 2 



Ej A 3 N(N + l)(N + 2) Xj, (83) 
is the NC electric field strength corresponding to NC Coulomb potential (the 
details will be published, see [IS]). 

We stress that in the NC case the electrostatic energy of electron, deter- 
mined by the trace in (1821) . is finite (no cut-off at short distance is needed). 
A straightforward calculation of the trace in (182]) gives the relation: 



mc 2 = | ^ A = | = A . (84) 

8 A 8 mc 1 

This Ao is fraction of the classical radius of electron ro = e 2 /mc 2 : Ao = 
1.06 x 10~ 15 m = 1.06 fm (the coincidence with the proton radius is purely 
accidental). 

The NC corrections to the if-atom energy levels given in this paper are 
of order (A /a ) = (9/64) a„ « 4 x 10~ n (here a ~ 1/137 is fine structure 
constant). Such tiny corrections to energy levels are beyond any experimen- 
tal evidence. Moreover, at Ao ~ 1 fm relativistic and QFT effects become 
essential. 

Our investigation indicates that the noncommutativity of the configura- 
tion space is fully consistent with the general QM axioms, at least for the 
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if-atom bound states. However, a more detailed analysis of the Coulomb 
problem in would be a desirable dealing, e.g, with the following aspects: 

• dyon problem (electron in the electric point charge and magnetic monopole 
field), Pauli if- atom (non-relativistic spin); 

• Coulomb problem in and its dynamical symmetry, QM supersym- 
metry and integrability of the Coulomb system. 

Besides non-relativistic if -atom, there are other systems that would be inter- 
esting to investigate within NC configuration space Rf , e.g. Dirac ii-atom 
(relativistic invariance?), non-Abelian monopoles, or spherical black- holes. 
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Appendix A 



Here we prove two formulas (1351) we need for the calculation of the NC 
Coulomb Hamiltonian. (We skip indices j and m here.) 
a. Let us begin with the first formula: 

[al[a a ,* jm ]] = \>[al K, J2 Iff : K : )] 

L — 4 m\ \ m2 ] . nil U2 ] - 

<J m ) 

= v E[a „M)!^4r ][a t, :K , <(-.)- 



mi! 77^2! Q ' 77J 712! 



z — ' mi! m2! 774! 772! 



(im) 



mi! 7712! Tlx! 772! 



+ A mJmJ ^' t 00 ' : U ^ 

^ (aj)^ (a 2 )^ t ag (-ag)^ 

+ A J > j j — [a a ,:K:][a J a , j — , (85) 

L — ' mi! m2! 77i! 772! 

(jm) 

where 71 = YlkLo c kQ k = SfcLo c kX k N k . Now we shall use the following 
commutation relations 

[4,:iV fc :] =-kal : N k ~ l : [4, : K :] = -Aa£ : : , 

[a a , : iV fc :] = jfe : A^" 1 : a Q [a a , : K :] = X : d e K : a a , (86) 

where <9 e denotes the derivatives with respect to £>: <9 e 7^ = ^ c fc f? fc l - 

It is easy to see that the second line in (|85|) vanish, and the the first and 
third line give the same contribution 



E (*^ ( _ Aj:W <l(^£l, (87) 

m\\ 777,2! nil n 2 l 

(J m ) 
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From (18T)]) the double commutator [a^ , [a a , : 7?. :]] follows directly, and this 
gives the value of the third line in f )85|) 

y- ^ l 2 _ A . 02^. +2A 1 \ 2 . 88 

' m\\m2\ tl\\ U2 ] - 

(jm) 

The last two equations yields the first formula in ( |35il . 
b. From equation (|27|) it follows easily 

N : N k : = : N k+1 : + A; : iV fc : £ : TZ : = : gTZ : + \: g d e TZ : . (89) 

This relation gives directly the second formula in (J35J) : 



(N+1)Y yUl) , l " 2 

' ^ m , I m.nl 



mi! m 2 ! 712! 
= J] ... [(iV + j + 1) :K:} ... 

= ^ ... : [{N + j + l)K + gd e TZ] : ... , (90) 

(jm) 

where we have replaced both untouched factors containing annihilation and 
creation operators by dots. 

Appendix B 

This part of the Appendix deals with how are normal and usual powers of 
g correlated, and how can this information be used to rewrite various kinds 
of power series. So the key relation of this section relates : g n : and g n . 

. „k . _ \k_N\_ _ I \\k( /\A . n -k . _ \ -fc N\ _ \-k 1 

■8 ■ ~ A W^hy. ~ y~ A ) -8 ■ ~ A jw+ky. ~ Wh 

(91) 
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(N) k stands for Pochhammer symbol. Since g = \N, there really are powers 
of g on the right-hand sides. 

Now to the above mentioned power series. Let us start with finding out how 
does the function : e@ e : modify when we rewrite : g n : in terms of g n in the 
corresponding Taylor series [fi denotes some arbitrary constant here): 



Considering the limit A — > the above equation corresponds to the knowm 
Euler's formula. A potential doubt arising from the colon marks on the left 
hand side ought to be dismissed due to the fact that : g n : — > g n if A — > 0. 
There is no way to distinguish between the normal and usual ordering in 
A = world. 

We can move to the more complex tasks now. In the course of many calcu- 
lations we need to handle expressions of the kind : g n e l3e : 



. p p e ■ - Si ■ n k ■ - V 



k=0 k\ (N-k)\ 



(92) 



(l + A/3) JV = (l + A/3)^ 



: g' l e 



n J3q . 



Ecx 
fc: 



oo /3 k :g' 
fc=0 fc 




E 



fc=0 fc! (AT-(n+fc))! 



= A 



(N-n)\ 



Nl 



E 



oo p k \ k {N-n)\ 
k=0 ~kT (N-n-k)\ 



= A 



AM 



(1 + /3A) 



(N—n)\ 
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And in the case of negative powers: 



e • l^k=0 k\ l^k=0 k\ (AT-(fe-n))! 



-n N\ Y^oo fi h \ k (N+n)\ 
(N+n)\ 2^k=0 k\ (JV+n-fe)! 



= (i + ?\r +n 

That is almost all that is needed to get rid of the normal powers in the 
solutions of the NCQM equations. Since hypergeometric functions are often 
written in the form of power series with coefficients expressed in terms of the 
Pochhammer symbols, it is useful to take notice of the following identities 
when it comes to handling factorials: 



(o) B =a(a + l)...(a + n-l) = ^f 



TV! N\ 



(N+n)l " (N+n)(N+n-l)...(N+l)N\ (N+l) n 



_ _ = N^N-n^n, = ( _ 1)B( _ ^ { _„ + + „ _ ^ 

= (-in-N) n 



Now we are going to supply the reader with some of the calculations which 
were skipped in the previous sections. To keep reasonable length of the for- 
mulas, we will replace the long expressions with a, c, j3, Q... substituting the 
"right arguments" instead of our abbreviations and completing the calcula- 
tion that way is straightforward... yet not so tempting. 
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The following derivation was left out in ( 15 ip : 



a; c; Qg) : 



: e 



Pa 



(q)m (Qg) m 
"i=0 (c) m ml 

^ n i \ oTN—m (o) T 



E' 



^m=0 



JV v^oo 



[i + m e 



m=0 



(c), 
QA 
l+A/3 

-QA 
l+A/3 



m . . 



AT! 



(c) m w! 



(JV— m)!m! 

i / l) m (-jV), 



(93) 



(a)m(-JV)m 1 

(c) m m! 



[l + A/3] iV F(a,-iV;c; ^) 



This is to be done and suitably used if one wishes to rewrite (]5Tj) in terms 
of the fundamental system (03]) and to get rid of the normal ordering only 
thereafter: 



: e^ip (a; c; Qg) : 



= : e Pe x 

v-oo (-l) m {a) m (a-c+l) m (Qg)- a - 
A Z^m=0 m! 



[1 + A/3] 



l+A/3 
AQ 



JV! 



X 



v^oo (a) m (a-c+l) 
Z_/m=0 m! 



(JV+a)! 

l+A/3' 



-AQ 



(JV+a)! 
(JV+m+a)! 



[1 + A/3] 



JV 



l+A/3 
AQ 



JV! 



(JV+a)! 



X 



oo (a) m (a— c+l) m 
m=0 (JV+a+l) m m! 



l+A/3 
-AQ 



[1 + A/3] 



JV 



l+A/3 
AQ 



JV! 



(JV+a)! 



X 



xF a, a — c+1, jV + a + 1. 



l+A/3 
-AQ 



(94) 



To complete (155]) . this was needed: 



„0O (_l)m + l( 2a - ) m+j + l/2 _ 



(2a) 



i+i/2 



E 



(2jr'+l+m)!m! 
oo (-l) m+1 (2aA) m (2j+l)!JV! 
m=0 (2jr'+l+m)!ra!(JV-m)! 



(2j+l)! 

(2a)^+ 1 / 2 v^oo (-JV) m (2aA)" 
(2j+l)! 2^m=0 (2j+2) m m! 

(2ap+ 1 / 2 



(2i+l)! 
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[-N, 2j + 2, 2a\) 



(95) 



The equation (1361) required: 



■ e -2 e /A -j-1/2 j .(VgT™) ■ = (~ 1 ) m+1 ( 2 ») m+J+1/2 ■ e" 2 ^™ ■ 



(2ay +1 / 2 ^00 (-l) m (l-2) iv (2aA) m (2j+l)!A^! 



(2j+l+m)!m! 

(2j+l)! ' ^m=0 (l-2)'™(2j+l+m)!m!(7V-m)! 
(2 Q y+ 1 / 2 (-l) JV + 1 v-00 (-N) m (-2a\) m 
(2j+l)\ Z^m=0 (2j+2) ra m! 

^^^"^ (-iV, 2j + 2, -2aA) 

(96) 



Appendix C 

In this section we are going to look on the scattering case more closely, 
namely the part dealing with decomposition of the solution fl67|) into two 
parts corresponding to the in- and out- going spherical wave. Firstly let us 
remind the form in which we got ( 1671) : 

R Ej = - F j + l-i-,-JV, 2j + 2; 2iAp ^ 

If we , for the sake of brevity, denote 

a = j + l-i2 , 6=-JV , c = 2j + 2 , z = 2Upg§f, (97) 
then can be written as 

R Ej = (1 _ z ) fo / 2 F (a, 6, c; z) (98) 

According to Kummer identities (see [16J), F(a,b, c;z) can be written as a 
linear combination of two other solutions of the hypergeometric equation 
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(I53"|) . namely 



(-z)- a F(a, o + l-c, + I-6; z" 1 ) 
and 

(z) a ~ c (l - ^) c - a - fc F(c - a, 1 - a, c + 1 - a - b; Z - 1 (z - 1)) (99) 

Take notice of the fact that if we decompose our radial solution ([51]) before 
handling the normal ordering, i.e. if we write the confluent hypergeometric 
function (|42p in terms of the fundamental system (|43p of confluent hypergeo- 
metric equation (T4"T|) . and deal with the normal ordering only thereafter (see 
Appendix B, eq. (194]) ). we will end up with ( 199]) again. So after using one of 
the numerous Kummer relations listed in [16J, we have 

R Ej = (1-zf 2 F(a, 6, c; z) = 

= r in(c-a) r(c)r(i-b) 
c r(o)r(c+l-a-6) 

x(^) a - c (l -^) c - a - fe / 2 F(c-a, 1-a, c+l-a-6; z- 1 (z-l)) + 

, „»7r(l-a) r(c)r(l-6) 
" r e r(c-a)r(o+l-6) 

x (-z)~ a (1 - z) h/2 F(a, a + l-c, a + 1 - 6; 2" 1 ) 

(100) 
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Substituting back (1971) leads to 

Re 3 = (^|) JV F( J + l-, f ,-iV ) 2 J +2;2^^|) 

p ( p-i\E \ 
\p+iXE J 

xF(j + l + zf, -j + z^ r+j + 2 + zf; -^ggf) 

^^IXe y x 
(101) 

Keep in mind that our endeavour is to compare the NCQM result with the 
QM one. It is therefore suitable to rewrite the prefactors in front of the 
hypergeometric functions as exponentials. This way we obtain, besides other 
factors, also logarithms of T-functions, which can be rewritten using the 
following formula (see [16]): 

In (T(z + a ))=(z + a-±\ ln(z) - z + \ ln(2vr) + £ ^^f n (102) 

\ / n=l ^ ' 

Here B n (a) is a Bernoulli polynomial. After certain rearrangements we finally 
acquire the result (1691) . 
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